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ABSTRACT 

We propose that the observed misalignment between extra-solar planets and their hot host stars can 
be explained by angular momentum transport within the host star. Observations have shown that this 
misalignment is preferentially around hot stars, which have convective cores and extended radiative 
envelopes. This situation is amenable to substantial angular momentum transport by internal gravity 
waves (IGW) generated at the convective-radiative interface. Here we present numerical simulations 
of this process and show that IGW can modulate the surface rotation of the star. With these two- 
dimensional simulations we show that IGW could explain the retrograde orbits observed in systems 
such as HAT-P-6 and HAT-P-7, however, extension to high obliquity objects will await future three- 
dimensional simulations. We note that these results also imply that individual massive stars should 
show temporal variations in their v sini measurements. 

Subject headings: internal gravity waves, angular momentum redistribution, extra-solar planets, hot 
j up iters 



1. INTRODUCTION 

Jupiter-mass planets with a few day periods (hereafter 
hot Jupiters) are found around 1-2% of solar type stars. 
At their present-day location ( < Q.Q5AU), in situ f or- 
mation through core accretion (Pollack et al. 1996) is 
challenging because al) their natal disks are relatively 
hot (> 10 3 K), a2) the local Keplerian velocity is > 10 2 
km s _1 (ie much larger than the surface escape speed of 
any known planets), and a3) the tidal perturbation of 
their host stars is strong. Consequently, bl) only a mod- 
est amount of refractory grains can avoid evaporation in 
this region, b2) high-speed collisions among planetesi- 
mals generally lead to their disruption, b3) the asymp- 
totic mass of protoplanetary embryos is limited, and b4) 
the gas accretion is quenched. 

These issues have led to a paradigm shift that hot 
Jupiters may have relocated over l arge distances from th e 
preferred sites of their assemblage Qlda fc Lin|2004||2010 ). 
The physical cause of planetary relocation remains a con- 
troversial issue. A widely adopted "disk-migration" sce- 
nario attributes the cause of hot Jupiters' orbital ev olu- 
tion to their tidal interaction with their natal disks (Lin 
et al. 1996[ ). Their orbital decay is stalled when they 
enter into a ce ntral cavity in duced by the host stars' 
magnetosphere ( |Koenigl||1991[ ) or when their natal disk 
is severely depleted. This disk-migration scenario natu- 
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rally accounts for the origin of multi ple planets which a re 
locked in mean motion resonances (Lee & Peale 2002) 

Distant gas giants may also venture close to their host 
stars if their eccentrici ty is n e ar uni ty through dynam 
ical excitation by the iKozai l (1 1962 1) resonance dWu 



Murray 2003; Nagasawa fc Ida 2011), planet-planet scat 



tering ([C hatter jee et al. 2008), or secular chaos (Wu & 
Lithwick' 201ip . Tidal dissipation within the scattered 



planets would circularize their orbits if their periastron 
distance can reach within a few times their Roche radii 



( |Rasio fc Ford 1996). However, mass loss during the 
extremely close encounters also leads to either escape 
( Guillochon et al.|2011 ) or res idual cores as short-period 
super Earths ( |Liu et al.||2Q12 ). This dynamical scenario 
predicts a pile up in the period distribution which does 
not seem to match t hat observed with a control sample 
flHoward et aT|2010| ). 

Recent observati ons o f the Rossiter-McLaughlin (RM) 
effect (Ohta et al. 2005[ ) indicate that the spin of many 
relatively massive and hot (> 6, 300K) main sequence 
stars is misal igned with their planets' angular m omen- 
tum vector flWinn et aLl|2010|]Schlaufman||2010| ). This 
misalign ment is an expected outcom e ol the dynamical 
scenario ( |Fabrycky fc Tremaine 2007), although its corre- 
lation with the stellar mass requires assumptions on the 
poorly know tidal dissipation mechanism inside main se- 
quence stars ( Winn et aO2011| |Lai|2012 ). 

The observed spin-orbit misalignment would pose a 
challenge to the disk-migration scenario for the origin 



of hot Jupiters (Triaud et al. 2010 Winn et al. 2011) 
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if the host stars retain the angular momentum they ac- 
creted from the protostellar disks, out of which the plan- 
ets formed. Within the observational uncertainty, this 
assumption is consistent w ith that observed in Kepler 30 
( Sanchis-Ojeda et al.|2012| ), but it cannot account for the 
7 degree tilt in the Sun's spin axis with respect to the 
angular momentum vector of the planetary orbits. 

Stars generally form in a turbulent cloud complex 
where fragments undergo chaotic encounters as they col- 
lapse to form stars through disk accretion. The spin ori- 
entation of protostellar disks may re-o rient e pisodi cally, 
especially during the disk formation. [Bate| ( 2Q1Q[ ) sug- 
gested that this process may have led to the observed 
misalignment between planets' orbital angular momen- 
tum vector and the spin axis of their host stars as well as 
differential rotation. But this possibility may depend on 
whether 1) there is adequate mass in the disk to form 
planets, 2) there is sufficient residual gas to re-orient 
planets' orbital plane during disk depletion, and 3) there 
is enough coupling between the stars and their surround- 
ing disks to realign the spin of their outer envelope. 

In this paper, we consider another possibility: that 
differential spin may naturally arise in relatively mas- 
sive and hot stars even in the limit that they started 
with uniform spin. We show that internal gravity waves 
(IGW) generated at the interface between the convective 
core and radiative envelope transport angular momen- 
tum outward. These waves dissipate in a shell near the 
surface of the star forcing the surface (including the pho- 
tosphere) to rotate differently than the bulk of the star 
and to vary in time. 

2. RESULTS 

2.1. Model Setup 

In order to study the angular momentum transport 
by waves in a massive star we solve the full set of hy- 
drodynamic equations in the anelastic approximation. 
These equations are solved in two-dimensional cylindrical 
coordinates representing an equatorial slice of the star. 
A complete description o f the equations and numerica l 
method can be found in Rog ers fc Glatzmaierj ( |2QQ5[ ). 
The reference state model is that ol a 3Mq star and is 
calculate d using th e Cambridge stellar evolutionary code 
STARS dEggleton|[l97T] ) f 7 ] In this model the convective 
core occupies the inner « T4% by radius and the radiative 
envelope extends to the stellar surface at « 1.6i?e[^| Over 
the simulated domain, the density varies from 40gm/cm 3 
in the core to 10~ 3 gm/cm 3 at the surface. Here we report 
mainly on the results of a model in which the convecting 
core is initially rotating 1.5 times faster than the radia- 
tive envelope and the average rotation is approximately 
ten times solar, designated MDR. An additional model 
with no initial rotation, designated MUR, will be briefly 
discussed. 

2.2. Angular Momentum Transport by Internal Gravity 

Waves 



1 Although we use a 3Mq model here we expect the results to 
be qualitatively representative of any star with core convection and 
an extended radiative envelope. 

2 We do not in this first attempt include a surface convection 
zone for numerical efficiency. 



In Figure 1 (a-e) we show the vorticity at five different 
times during the simulation for MDR. We can immedi- 
ately see that the small convecting region efficiently gen- 
erates IGW0 which propagate into the radiative envelope 
with amplitudes similar to that of convection (wave ve- 
locity amplitudes at generation are « 10 3 — 10 5 cm/s). 
The wave generation process is extremely time and lo- 
cation dependent. In Figure la one can clearly see that 
waves are generated where an upwelling plume hits the 
convective-radiative interface. At one time waves are 
generated at one main location, at another a different 
location and at another at multiple locations (see Fig- 
ures lb and lc as well). Therefore, the wave pattern 
seen in the radiative envelope is some combination of 
waves generated at various times and locations. One can 
also see in Figure la that the wave patterns are predom- 
inantly horizontal or concentric, reflecting the fact that 
these waves propagate mostly horizontally because of the 
smallness of cu/N where uj is the wave frequency and N 
is the Brunt- Vaisala frequency. The wave propagation is 
therefore, a spiral propagation outward. 

Once generated, the propagation and dissipation of 
the waves depends sensitively on the properties of the 
waves. Low frequency waves are generated with higher 
amplitudes than high frequency waves (the generation 

spectrum we measure is A(uj) oc e ( blUJ + b2UJ )). How- 
ever, on their propagation outward low frequency waves 
are dissipated by radiative diffusion more strongly than 
their high frequency counterparts. The radiative diffu- 
sion acts nearly as an equalizer, so that while high and 
low frequency waves are generated with different ampli- 
tudes, once they reach the surface there is less disper- 
sion. The waves that reach the surface with the highest 
amplitude are those waves with high enough frequency 
to escape substantial radiative diffusion but low enough 
frequency to be generated efficiently, these tend to be 
waves with frequencies « lO/aHz and large scales (hor- 
izontal wavenumber, k, 1-3). It is these mid- frequency 
range waves that contribute the most to the angular mo- 
mentum transport at the stellar surface. Waves with 
this combination of frequency and wavenumber reach 
the surface with velocities approaching « 10 5 cm/s, have 
Froud^] numbers close to 1 and are therefore, nonlinear. 
Some evidence of wave breaking, rippled surfaces and 
large amplitudes, can be seen in figure lb. While we see 
some evidence of wave breaking, the actual process by 
which waves initially dissipate in the outer layers is not 
fully understood. Several physical processes could con- 
tribute including: 1) nonlinear wave breaking 2) wave- 
wave interaction as waves reflect off the top boundary 
or 3) wave- wave interactions as waves internally reflect 
when propagating into a region of rapidly varying N. 

Because of the initially imposed differential rotation 
at the convective-radiative interface, prograde waves are 
doppler-shifted to higher frequencies, while retrograde 
waves are shifted to lower frequencies. Therefore, ra- 
diative diffusion, being frequency dependent (the radia- 
tive damping length is proportional to a; 4 ), dissipates 

3 We use the term "wave" here loosely to mean "disturbance". 

4 The Froude number is a general measure of wave linearity and 
is defined as w'/LN, where w' is a velocity amplitude, L is the 
wavelength and N is the Brunt- Vaisala frequency. Numbers close 
to unity indicate nonlinearity. 
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retrograde waves slightly more than prograde waves. 
Consequently, upon arrival at the surface, retrograde 
waves have slightly lower amplitude than their prograde 
counter-parts and therefore, contribute slightly less to 
the angular momentum evolution. This means that, al- 
though both prograde and retrograde waves reach the 
surface and dissipate, over time, prograde waves deposit 
more angular momentum and the region accelerates. 

Eventually the shear flow reaches an amplitude such 
that it represents a critic al layer to many of the d omi- 
nant phase speed waves ( |Booker fc Brether ton 1967]). A 
critical layer is defined as the position where the local 
angular velocity is equal to the horizontal phase speed of 
the wave, c p h- In the ray-tracing description of IGW 
the vertical wavenumber goes to infinity at a critical 
layer, therefore, the time to approach the critical layer 
goes to innnityjj This means that any small dissipation 
leads to complete absorption of the wave at the criti- 
cal layer. For example, at the position labeled (1) in 
Figure 2b, the angular velocity is « 10 _5 Hz, therefore, 
a wave with frequency 10/iHz and horizontal wavenum- 
ber one, has a horizontal phase speed 10 _5 Hz and ex- 
periences critical layer absorption. Indeed, any wave in 
which c p h = uo/k > 10 _5 Hz will encounter the critical 
layer and deposit most, if not all, of its angular momen- 
tum locally. This causes the rapid increase in the angular 
velocity seen between times labeled (1) and (2) in Figure 
2b. 

To quantify the efficiency of angular momentum trans- 
port by IGW we turn to the evolution equation for the 
mean angular velocity: 



2a), si milar to what occurs in the Quasi-Biennial oscil- 
lation (Bald win et al.pOOl ) in our own atmospher e and 



to the Plumb fc M cEwan laboratory experiment ( Plumb 



dU _ 1 dr(pWH) | ^ (d 2 U 
dt pr dr \ dr 2 



r p J dr 



the first term on the right hand side (rhs) of (1) repre- 
sents the angular momentum transport by IGW or the 
wave momentum flux (the second term on the rhs of (1) 
represents viscous dissipation, near the surface the vis- 
cous dissipation is always at least two orders of magni- 
tude smaller than the momentum flux by waves). This 
value is ~ 10cm / s 2 at the base of the radiative envelope 
and around lcm/ s 2 at the surface, although there are 
order 10-100 variations in time. This can lead to rapid 
time variation of the mean flow. The timescale for the 
prograde layer to grow to its peak value is extremely 
rapid. Within about 9 rotation periods, or approx- 
imately lO -10 ^^, where tes is the Eddington- Sweet 
time, the angular velocity at the surface has increased 
by a factor of 10. Stated another way, in the rotating 
frame of reference the linear velocity at the surface has 
increased from zero to « 10 6 cm / s in 9 rotation periods, 
consistent with the values of momentum flux given above. 

As time progresses, the prograde layer migrates as con- 
tinuously generated waves steepen and dissipate within 
the critical layer. This causes the prograde layer to 
propagate toward the convection zone in time (Figure 

5 The critical level is defined as the place where the horizontal 
phase speed of the wave is equivalent to the mean flow speed. As 
the wave is incident on the critical level there is a region below 
the critical level where substantial wave dissipation occurs, which 
is referred to as the citical layer. We have not necessarily adhered 
to these subtle distinctions in our physical description. 



fc McEwan| [l978). Eventually the entire radiative enve- 
)e is spun up. However, beca 



lope is spun up. However, because of angular momentum 
conservation, the convection zone slows down (see Fig- 
ure 2a). This produces a differential rotation profile that 
favors retrograde waves at the surface, so the process re- 
verses and eventually, the surface will spin retrograde. 
We have not been able to run this model long enough for 
the reversal to occur, however, one can see the slow decay 
of the prograde shear in Figure 2b as well as retrograde 
pockets in Figure le. At the current rate of deceleration 
(which is substantially slower than the rate of accelera- 
tion since no retrograde critical layer has yet formed) the 
prograde radiative envelope would take approximately 
a few hundred rotation periods to reverse, or lO -7 ^^, 
which astronomically speaking, is still an extremely short 
timescale. The growth and reversal timescales are in- 
verse ly proportional to the wave momentum flux ( |Plumb| 
1977). The reversal timescale is longer than the initial 



growth time because the convection zone contains most 
of the mass, therefore, when it slows down the differen- 
tial rotation between it and the surrounding radiation 
zone is weaker than the initially imposed differential ro- 
tation. Hence, the amplitude variation between prograde 
and retrograde waves at the surface is reduced and the 
net momentum flux is lower. 

In Figure l(f-j) we show the vorticity for the model 
MUR, in which there is no initial rotation. Here we 
clearly see dynamics similar to that seen in Figure (a-e) 
but the initial surface shear that develops in this model is 
retrograde. In a model with uniform initial rotation the 
sense of the initial surface shear layer is random as there 
is no preferred direction. Otherwise, this model evolves 
similarly to the one previously described with the devel- 
opment of a critical layer, rapid momentum transport 
and momentum fluxes similar to those quoted above for 
MDR. Although we have only presented these two models 
we note that this is a very robust mechanism as a sur- 
face shear flow has developed in all models we have run, 
whether initially uniformly or differentially rotating. The 
amplitude of the shear and its growth rate depend on the 
model details but the time dependent shear development 
is rigorous. The two models presented demonstrate the 
large variation in amplitude and sign that surface flows 
on hot stars can experience. 

3. DISCUSSION 

In the above models we have shown that IGW can 
efficiently transport angular momentum in stars with 
convective cores and extended radiative envelopes. This 
transport causes the angular velocity profile to vary sub- 
stantially in radius and time. In the model presented 
the angular velocity at the surface changed by a factor 
of 10 in a timescale minuscule compared to other rele- 
vant timescales. In all cases run the angular velocity 
changes in time. This has significant consequences for 
observations of massive star rotation, as these models 
indicate that the rotation of the surface layers is not nec- 
essarily representative of the rotation of the entire star. 
They also indicate that the surface rotation of the star is 
time dependent. This is a robust process, whic h happens 
in our own atmosphere ( Baldwin et al.pOOT ), has been 
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demonstrated experimentally ( Plumb fc McEwan| [l978), 



and shown here numerically. It likely occurs in any star 
with a convective core and extended radiative envelope 
and has many interesting consequences. 

This mechanism provides a natural explanation for 
why mis-aligned planetary systems tend to be around 
hot stars. In hot stars which have convective cores and 
radiative envelopes IGW transport sufficient angular mo- 
mentum to cause the surface of the star to rotate differ- 
ently than it did when the planetary system was formed. 
Although our 2D models only strictly explain retrograde 
systems we view this work as a proof of concept and 
application to other high-obliquity objects will await fu- 
ture 3D simulations^] These simulations imply, at the 
very least, that the equatorial angular velocity at the sur- 
face could be time-dependent as the IGW-driven shear 
flow varies in time. This physical effect would lead to 
variations in measured RM effect. Multiple observation s 
of the exoplanet sy stem XO-3 bylHebrard et al. (2008), 
Winn et ah] ( |2009| and |Hirano et al.| ( |2Ulip reveaTan 
unexplained variation in spin-misalignment , which could 
be interpreted as a 30% variation in equatorial velocity 
over a year. Multiple observations of the RM effect in the 
same hot-star system over time could provide an impor- 
tant test. Similarly, inclination measurements of massive 
stars over time should show variations. 

One way to distinguish this process, in which the 
mis-alignment is rooted in angular momentum varia- 
tions in the star alone, from one in which the plan- 
ets themselves are mis-aligned, would be the observa- 
tion of multiple planets in the same system with simi- 
lar mis- alignments. This would be difficult to produce 
in a scattering-circularization scenario. Recent obser- 



vations of the mu lti-planet system Kepler-30 (Sanchis- 
Ojeda et al.|20 12) around a solar type star does not pro- 
vide a test of our model. Observations of multi-planet 
systems around a hot star would provide a crucial test. 

This wave-mean flow interaction could have several 
other interesting implications. The surface angular ve- 
locity achieved in this model is unstable to the shear 
instability, likely resulting in turbulence and substantial 
mixing of species. The angular velocities achieved here 
are about 2%Q C and higher angular velocities are likely 
possible. These amplitudes, combined with the short 
timescale of angular velocity variation could explain the 
observed BE class of stars. 

How this mechanism operates in three-dimensions (3D) 
is not clear. In our own atmosphere IGW driven shear 
flows are limited to equatorial regions because the IGW 
forcin g is offset by the C oriolis force at higher lati- 
tudes ( [Baldwin et al.pOOl ). Although the analysis may 
be slightly altered in massive stars (as hydrostatic and 
geostrophic balance do not apply), it is highly likely these 
stars have latitudinal differential rotation, which should 
be accounted for when modeling the RM effect. More 
detailed analysis of latitudinal and temporal variation in 
surface angular velocity await more sophisticated simu- 
lations. 

We suggest several observational tests of this scenario. 
RM observations should be analyzed assuming some lat- 
itudinal differential rotation. Repeated observations of 
the RM effect associated with the same hot star sys- 
tem would be a critical test and finally, observation of a 
multiple-planet system around a hot star in which all of 
the planets had the same mis-alignment would provide 
strong support for this model. 
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Fig. 1. — (a-e) Time snapshots of vorticity at five different times during the simulation MDR, corresponding to (a) 4.8 x 10 5 s, (b) 
1.76 x 10 6 s, c) 4.6 x 10 6 s, (d) 2.2 x 10 7 s and (e) 4.4 x 10 7 s. (f-j) Time snapshots of vorticity at five different times during the simulation 
MUR, corresponding to (f) 2.57 x 10 6 s, (g) 2.66 x 10 6 s, (h) 5.08 x 10 6 s, (i) 1.15 x 10 7 s and (j)1.23 x 10 7 s. White represents positive 
vorticity, while black represents negative vorticity. 

in Solar physics under award number 0457631. D.N.C. Lin was supported by NASA (NNX08AM84G), NSF (AST- 
0908807) and University of California Lab Fee. Computing resources were provided by NAS at NASA Ames. 



Rogers et al. 




Fig. 2. — (a) Angular velocity as a function of time (x) and radius (y) for model MDR. Red represents prograde flow while blue represents 
retrograde flow. The blue horizontal line represents the top of the convection zone. We see a prograde layer form initially at the stellar 
surface and then propagate toward the convection zone in time. Initially the convection zone and overshoot layer show both prograde and 
retrograde flow. However, as the radiative envelope is spun up, by angular momentum conservation the convection zone must slow, as is 
seen, (b) Line plot of angular velocity in time for radii (moving from top to bottom) 0.98, 0.86, 0.45 and 0.17 of the stellar radius. Angular 
momentum deposition occurs by critical level absorption between time labeled (1) and (2). The angular velocity has increased by a factor 
of 10 in just 9 rotation periods. Subsequent decay represents the gradual reversal process as retrograde waves now reach the surface with 
slightly higher amplitude. 



